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Quasi-static systems are an important concept in thermodynamics: they are dynamic but close 
enough to equilibrium that many properties of equilibrium systems still hold. Slowly evolving 
horizons are the corresponding concept for quasilocally defined black holes: they are "nearly isolated" 
future outer trapping horizons. This article reviews the definition and properties of these objects 
including both their mechanics and the role that they play in the fluid-gravity correspondence. It 
also introduces a new property: there is an event horizon candidate in close proximity to any slowly 
evolving horizon. 
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When trying to understand a physical system, one often starts with its equilibrium states. For example, introduc- 
tions to mechanics generally begin with statics while the zeroth law of thermodynamics defines a system in thermal 
P-h equilibrium. After the time-independent physics is understood, one turns to dynamics. However this generally entails 
a significant leap in complexity: thermodynamics is over 150 years old but non-equilibrium thermodynamics is still an 
^vq active area of research. Thus as an intermediate step in the approach to full dynamics it is often useful to start with 
the near-equilibrium regime where many of the lessons from equilibrium physics still apply. For thermodynamics, this 
1 ^ 1 is the quasi-static regime where a system smoothly (and generally slowly) transitions between equilibrium states. 

The study of black holes follows the usual pattern. In four-dimensional asymptotically flat space-times, the equi- 
librium states are identified with the Kerr-Newman family of exact solutions. Well-known theorems identify them as 
the unique stationary, axisymmetric and asymptotically flat vacuum black hole solutions. As such they have been 
i— — i intensively studied since their discovery almost 50 years ago. Much of what we know about black holes comes from a 
study of these equilibrium solutions and their near-equilbrium perturbations (see, for example, [T]). 

That said, as has been repeatedly emphasized elsewhere in this volume, this family of solutions is not the end of 
the story. Real black holes don't sit alone in an otherwise empty universe and so there is great interest in developing 
qq a more localized theory which characterizes black holes by their physical and geometric properties rather than as any 
[~^. particular exact solution. Hence the closely related programmes studying apparent, trapping, isolated and dynamical 
horizons (review articles include [2 6 ). In this context, isolated trapping horizons are the equilibrium states of interest 
and dynamical trapping horizons are the corresponding time-dependent states. 

In this chapter we will review the geometric characterization and physical properties of near-equilibirum black hole 
£\J horizons. These are the "almost-isolated" slowly evolving horizons first defined by Booth and Fairhurst [7J [S]. In 
section Hlwe consider the geometrical background common to all types of horizons. With this material in hand, 
section ffl] reviews the various definitions of horizons, focussing on slowly evolving horizons. Section |IV| discusses 
•rH two of the key properties of these horizons: they obey dynamical laws of black hole mechanics and also are always 
accompanied by a candidate event horizon. Section [V] considers a couple of examples of the slowly evolving horizons: 
a simple Vaidya spacetime and some much more interesting five-dimensional black brane spacetimes that manifest 
the fluid-gravity duality. Finally section [VT| closes this chapter with a summary of our main points. 



II. BACKGROUND 



We begin with a brief review of the geometry of two- and three-dimensional surfaces in four-dimensional spacetime. 
Much of this material is repeated in other chapters of this volume, however this section will serve to establish notation 
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and perhaps present things with a slightly different perspective. For those who would like more details, this is 
essentially an abbreviated version of the extended discussion in [8]. 



A. Intrinsic and extrinsic horizon geometry 

Let (M, <7a&,V a ) be a four-dimensional, time-orientable spacetime and (H,q ab , D a ) be a three-dimensional hyper- 
surface which can be foliated into two-dimensional spacelike surfaces (S v ,q ab ,d a ). The normal space to each of the 
S v can be spanned by a pair of future-oriented null vectors £ a (outward-pointing) and n a (inward-pointing). The 
direction of these vectors is fixed, but their scaling isn't. One degree of freedom is usually removed by requiring that 
£ ■ n = — 1 which leaves a rescaling freedom: 

£^re s £ and n^e~ f n, (1) 
for an arbitrary function /. Independent of the choice of scaling, we can write the two- metric as 

q~ab = 9ab + £ a n b + n a £ b (2) 

and the four-dimensional volume element e as 

e = ^AnAe (3) 

where £ and n are the one-form versions of the corresponding vectors and e is the induced volume element on the S v . 

The intrinsic geometry of the S v is determined by the induced metric q ab but the extrinsic geometry is defined by 
how the null normals vary along the surface. We have the extrinsic curvature analogues: 

= FMVJU and *>> = q c a q d b V c n d (4) 
as well as the connection on the normal bundles: 

u a = -q b a n c V b £ c • (5) 
It is useful to decompose the extrinsic curvatures into their trace and trace-free parts: 

k ah = 2 ( ( )1ab + <J ab and k ab > = -6 (n) q ab + a ab ' , (6) 

where the traces are called expansions and the trace-free parts are shears. These all depend the scaling of the null 
vectors and under the rescalings defined by ([lj: 

e?k$ , e-fk { 2 and u a -> Q a + d a f . (7) 

Returning to H, one can further reduce the freedom allowed to the null vectors by tying their scaling to the foliation. 
Given a coordinate labelling v of the foliating two-surfaces, a unique future-oriented vector field V a on H satisfies the 
following conditions: 1) it is normal to the S v , 2) it is tangent to H and 3) it satisfies Cyv = 1. This is the evolution 
vector field for that foliation labelling: it evolves leaves of the foliation into each other. In coordinate language 
V = d/dv for any parameterization in which V° Lie-drags the spatial-coordinates between S v . Independently of 
whether or not such a system is constructed, the scaling of the null vectors can be fixed by the requirement that 

V a = t - Cn a , (8) 

for some function C which we call the expansion parameter. Note that if C > 0, H is spacelike, while if C < it is 
timelike and if C — it is null. We will mainly be interested in the cases where C > 0. 

Given this construction, the scaling freedom of the null vectors is restricted to the freedom to reparameterize the 
foliation labelling. For an alternative labelling v — v(v) we have 

[di] a = -7^[dv} a => V a = a(v)V a (9) 
a(v) 

where a(v) = || is constant over each S v . Then 



T = at , h a = -n a and C = a 2 C . (10) 
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Under this restricted class of rescalings with a constant over each individual S v , Cu a is invariant. 

The intrinsic and extrinsic geometry of the full H can be expressed mainly in terms of the two-surface quantities. 
First, the induced metric is 

q ab = 2C[dv] a [dv]b + q a b , (11) 
which explicitly demonstrates how the sign of C determines the signature of the metric. Next, defining 

r a = 4 + Cn a (12) 
as a future-oriented normal to the surface, the associated extrinsic curvature is 

K$ = q c a qlV c T d = (2Ck v - C V C) [dv] a [dv] b + 2C {Q a [dv] b + [dv] a Q b ) + (k ( 2 + CA^ } ) , (13) 

where Ky = — V a n b V a i h . Note that we have not unit-normalized r a so as to allow for all values of C. If H is spacelike, 
the usual extrinsic curvature define in terms of the timelike unit normal f a is 

K [f) - 1 K (t] (141 

B. Evolutions and deformations of two-surfaces 

Next, consider how these geometric quantities change if the surfaces are deformed. We start with the geometry of 
H for which the formalism will probably be more familiar as it identical to that used in the (3 + l)-decomposition 
of general relativity (see, for example, [HE]). Restricting our attention to normal deformations and a spacelike H, 
potential evolutions are generated by vector fields of the form T a — Nf a , where N is a function. Then 

5 T qab = 2NK$ (15) 

while 

5 T K%> = D a D b N - N [R ab + kWk$ - 2K^ C K% + q a c q b d K cd ) , (16) 

where as specified earlier D a is the induced covariant derivative on the surface, R ab is the (intrinsic three-dimensional) 
Ricci tensor and lZ ab is the Ricci tensor for the full spacetime. In (3 + 1) -general relativity these are the time-evolution 
equations for evolutions with lapse function TV but vanishing shift vector. 

Similar equations can be developed for deformations of the individual two-surfaces. It is probably easiest to visualize 
this with the help of a coordinate system. Relative to a full coordinate system {x a } on some region of M, we can 
parameterize any S v by coordinates 6 : so that S v is embedded in M by the relation x a = X a {9 A ) for four functions 
X a . Again we restrict our attention to normal deformations, noting that any normal vector field X a (6) defined over 
S v can be written as 

for some functions A and B. Then 

X a {6 A ) ^ X a (6 A ) + eX a {0 A ) (18) 

defines a new surface S' v by deforming S v a coordinate distance e in the direction X a . Further it identifies points on 
the two-surfaces (essentially by Lie-dragging coordinates between surfaces) as shown in Figure [I] Then, for example, 
the deformation of the two-metric q ab in the direction X a is defined by 

x ~ v qAB\ x+eX - qabIx 

OxqAB = hm , (19) 

where we have written the components of the metric in terms of the surface coordinates. 

As for the (3 + l)-decomposition of general relativity, computationally this operation amounts to taking Lie deriva- 
tives with some extra conditions imposed to ensure that the geometric quantities are calculated on the correct surface 



curve of constant 8 A 



FIG. 1: This diagram can be interpreted in two ways. First letting X a — V a it may be thought of as a full horizon H with 
the S v and S' v as foliation surfaces. Alternatively it depicts a more general deformation of an S v into S' v in which case X a is 
the deforming vector field. Of course the dual interpretation arises because the evolution of the horizon is just a special case 
of more general deformations. 



and that all quantities are appropriately pulled-back into those surfaces (full details may be found in [5]). Then it is 
straightforward to show that 



5 e q ab = 2k ( J and S n q ab = 2k { ^ 



(20) 



while 



See = ' 



->{£) 



and S„e = et 



(n) 



(21) 



Now it is clear why the traces of the extrinsic curvatures are called expansions: they tell us how the area element 
changes when the surface is evolved in the given direction. The shears are the part of the evolution that deforms the 
surface (but does not change its area). 

Normal deformations are linear for the first derivatives of the metric. For example: 



S x e = AS e e - B6 n e = e(A6 w - B0 [n) ) 



(22) 



for arbitrary functions A and B. However, as in the three-dimensional case, things become significantly more com- 
plicated for second derivatives. For example 



- d 2 B + 2u a d a B 



-B 



\wf - d a u a - R/2 + G ab £ a n» - 6 {l) 6 {n) 



A 



G ah ti b + (l/2)0fo 



(23) 



and 



S x 9 (n) = -K X e {n) +d 2 A + 2Co a d a A 



-A 



l,~.l|2 



+ d a u a - R/2 + G ab n a £ b - (i) (n) 



B 



k Cn) f + G ab n V + (1/2)0 



(n) 



(24) 



where k x = — X a n b W a £ , R is the two-dimensional Ricci scalar for S v , \\uj\\ 2 = w°w a , ||c(£)|| 2 = vfe^ab an d \\ a (n)\ 
a ?n) a ab ■ ^ 0r * ne connection one- form: 



5 X Cb a = -(A6 W - B0( n ))Cb a + d a n v - d b \ Aa^ a + Ba^ a J - q b G bc i 
+ 7jd a (A0(e) + B9( n) ) — 9^d a A — 6^ n )d a B . 



(25) 
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where, as before, r c = A£ c + Bn c is normal to X c . All of these second derivatives depend on how A and B vary over 
the S v . 

Probably the most intuitive application of these deformations is in understanding the time-evolution of a horizon. 
Then X a is identified with the V a evolution vector field and, for example, 

S V 9 {1] = k v 9 {£) - d 2 C + 2£o a d a C (26) 

-C [|H| 2 - d a u a - R/2 + G ab tn h - W (n) ] - [h W f + G ab £ a £ b + {l/2)6 2 {i) . 

In particular for a null surface V a = £ a and we recover the Raychaudhuri equation 

8 t e {l] = Kl 9 w - [\\a (e) || 2 + G ab £ a £ b + (1/2)0 2 J , (27) 

which is key to many of the calculations involving isolated [5] or event [TJ |5] horizons. 

III. GEOMETRIC HORIZONS 

We are now ready to define slowly evolving horizons. In this section we will use geometric arguments, leaving 
physical properties for the next section. We begin with a review of the definitions of isolated and dynamical trapping 
horizons. 

A. General quasilocal horizons 

In a four-dimensional spacetime (M, g ab , V a ) a future outer trapping horizon (or FOTH) is a three-surface H that 
is foliated by spacelike two-surfaces (S v ,q a b,d a ) such that on each surface: i) 9^ = 0, ii) 0(„) < and hi) there is a 
positive function f3 such that Sp n 9^ < [10 . These conditions are intended to (locally) mimic those used to define 
apparent horizons [IT]: each slice of a FOTH is marginally outer trapped (9^ = 0) but the other two conditions 
guarantee that it is possible to deform the S v inwards so that they become fully trapped. 

That the inward expansion be negative is straightforward to check and independent of the scaling of the null vectors, 
however the last condition is not quite so simple. From the general deformation equation for 9^ we see that with 
9 W =0: 

Sp n 6(i) = d 2 p - 2u a d a p - (R/2 + ||w|| 2 - d a Cb a + G ah £ a n b ) (3 . (28) 

Now, for a given scaling of the null vectors, deciding whether condition iii) holds will amount to a study of the possible 
behaviours of this second order elliptical differential operator. Even for fairly simple cases, a non-trivial f3 can be 
required. For example, rapidly rotating Kerr solutions and the standard scaling of the null vectors obtained from 
Boyer-Lindquist co-ordinates has 6 n 9^) it (appendix C of 8J). However, the horizons certainly are FOTHs and this 
can be demonstrated with a suitable choice of /?. Similarly for general FOTHs away from spherical symmetry, it will 
often be the case that the scaling determined by fixing V a = £ a — Cn a , is not a scaling for which 5 n 9^ < 0. Again a 
non-trivial (3 is required. 



Since 9^ — on a FOTH we much also have = 0. Then applying a maximum principle argument to (23) 

along with the fact that 6p n 6m < (for some /3) one can show that the null energy condition implies that C > 
which in turn means that 

C v e = -C9 (n) e > 0, (29) 

since < 0. Thus, there is a second law of FOTH mechanics: their area is non-decreasing. 

Next, we summarize isolated horizons (see [5] in this volume for a more information or [T^J for the full details). A 
three-dimensional submanifold in that same spacetime is a non-expanding horizon if: i) it is null and topologically 
S x IR for some closed two-manifold S, ii) 9^j — and iii) —T ab £b is future directed and causal. As usual £ a is a future 
oriented, outward-pointing normal; since the horizon is null, in this case no foliation is needed for its construction. 
That said a foliation is certainly no hindrance to a three-surface being a non-expanding horizon and any null FOTH 
satisfying the null energy condition will certainly be a non-expanding horizon. 

Non-expanding horizons are the simplest objects in the isolated horizon family. Any non-expanding horizon can be 
turned into a weakly isolated horizon if the scaling of the null vectors is chosen so that C gw a = where 

w a = -n b V a £ b = -Kgn a + u) a , (30) 
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and the arrow indicates a pull-back into the cotangent bundle of the non-expanding horizon. With this scaling, zeroth 
and first laws of isolated horizon mechanics may be established. Finally, if there exists a scaling of the null vectors for 
which the entire extrinsic geometry of the horizon is invariant in time then the weakly isolated horizon is an isolated 
horizon. Thus in that case not only do Liki and CiHj a vanish, but so do CiOt n ) and £i<j^ . All Killing horizons are 
isolated in this sense. 

Our final object is a dynamical horizon [T3], A three-dimensional sub-manifold of a spacetime (M, g ab ) is a dynamical 
horizon if: i) it is spacelike and ii) can be foliated by spacelike two-surfaces such that the null normals to those surfaces 
satisfy Q(j\ = and 6r n \ < 0. Perhaps the key property of dynamical horizons is that they always expand: the spacelike 
assumption means that C > and so by the strictly greater-than version of (29 1, dynamical horizons always expand. 
Not surprisingly, if the null energy condition holds and 6i9m ^ at least somewhere on each S v , then a FOTH is 
a dynamical horizon. Another important property of dynamical horizons is that their foliation into S v surfaces with 
9i{\ =0 is unique |14j . This is very convenient as we do not need to worry about whether or not geometric properties 
are foliation dependent: they are, but since the foliation is unique this is fine! 

Note that this foliation rigidity contrasts with the corresponding situation for isolated horizons where foliations can 
be freely deformed. However, the situation is very different if we consider the rigidity of the full three-dimensional H. 
Marginally outer trapped surfaces in an isolated horizon can only be deformed along the horizon itself. By contrast 
for a dynamical FOTH, the 0i(\ = surfaces can only be deformed out of H (this is essentially equivalent to how 
the exact location of a dynamical apparent horizon depends on the foliation of the full spacetime). For an extended 
discussion of this point see [5] . 

We will take FOTHs as our basic objects and classify them by hybridizing the naming systems as was done in [5]. 
Thus a FOTH that also satisfies one of these other sets of properties will be referred to as non-expanding, (weakly) 
isolated, or dynamical as appropriate. There are also more exotic forms of FOTHs associated with apparent horizon 
"jumps" (see for example dSJUl]), inner horizons or white holes [TU]. However, they can be left aside for the purposes 
of this article. The horizons that we will be interested in will always be dynamical or weakly isolated FOTHs. 



B. Near-equilibrium quasilocal horizons 

Intuitively, a near-equilibrium black hole boundary should be "almost" isolated. Thus, since isolated horizons are 
non-expanding and null, a near-equilbrium horizon should expand slowly and be almost null. The complication here 
is quantifying "expand slowly" and "almost null" : expanding dynamical horizons are, by their very nature, spacelike 
so it is not immediately clear what "slowly" means and further it is certainly not obvious how a spacelike surface can 
be "almost null". 

In analogy with isolated horizons, we break the definition up into two parts: slowly expanding horizons which focus 
on restrictions to the intrinsic geometry and slowly evolving horizons which add constraints to the extrinsic geometry. 



1. Slowly expanding horizons 



Definition: Let AH be a section of a future outer trapping horizon foliated by two-surfaces S v so that AH = {U V S V : 
vi < v < U2}- Further let V a be an evolution vector field that generates the foliation so that Cyv = a(v) for some 
positive function v, and scale the null vectors so that V a = £ a — Cn a . Finally let Rh be the characteristic length scale 
for the problem. Then AH is a slowly expanding horizon if the dominant energy condition holds and 

1. \R\,Cb a Cb a , \d a u a \ and T ab i a n b < l/R%. 

2. Two-surface derivatives of horizon fields are at most of the same magnitude as the (maximum) of the original 
fields. For example, ||d a C|| < C max / Rh, where C max is largest absolute value attained by C on S v . 

3. e <C 1 where e 2 / R 2 H — Maximum 

If the a(v) is chosen so that C « e 2 then this scaling of the null normals is said to be compatible with the evolution 
parameter e. 

This definition requires some discussion and we analyze it point-by-point. Starting with the preamble, we specify 
a section AH rather than a full horizon so that the definition will cover sections of horizons that may be slowly 
expanding only for finite periods of time. For standard four-dimensional black holes the characteristic length scale 
Rh will be set by the areal radius Rh = \/A/{4<k). 



C (V")|| 2 + T ab n a n b + 9f n) /2) . 
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Next, following from the discussion back in section II A, choosing Cyv = a(v) for some positive function a(v) 



ensures that the flow generated by V a evolves leaves of the horizon into each other and so can be thought of as a 
coordinate vector field compatible with the foliation. The free function is included so that the scaling is not tied to a 
particular choice of foliation labelling; including it is equivalent to allowing for relabellings of the S v . 

Turning to the main clauses, the notation X < Y means that X < k Y for some constant k of order one. Then 
our first condition that 



I and \G a b 



< 



1 

R H 



(31) 



is essentially a restriction on the allowed geometries of the horizons. The first part ensures that the curvature is not 
too extreme, the second two assumptions place some restriction on the extrinsic curvature while the last bounds the 
Einstein tensor (or equivalently stress-energy). In particular it was shown in [5] that these restrictions hold for all 
members of the Kerr family of horizons. An immediate consequence of these assumptions is that by (28) 



Hi) 



< 



1 



(32) 



is similarly bounded. Note the use of the absolute value sign. Even though we must have 5p n 6m < for (3 scaling, 
in general this will not be true for (5 = 1. 

The second condition restricts the size of the derivatives of the horizon fields so that they are at most commensurate 

(n) 

with their maximum values. For tensor quantities such as a J the magnitudes are taken to be judged relative to an 
orthornormal frame. For example 



>) 

J bc 



< 



1 



R 



H 



T (n) 
J bc 



(33) 



indicates that this inequality should hold for all components relative to any orthonormal frame on S v . 

We can then consider the last condition. After the set-up of the earlier conditions, this is the one that captures the 



essence of what it means to be slowly evolving. Beginning with the evolution of the area- form e on S v , by ( 22 ) we 
have 



C ye 



-CQ( n )€, 



(34) 



which is certainly dependent on the scaling of the null vectors. However this dependence may be easily isolated by 
rewriting 



C v e = -C0,„<Z = IIVII - 



(35) 



so that the effects of the rescaling freedom are restricted to pre- factor ||V|| = V2C. The term in parentheses (which 
is required to be small by the third condition) then provides an invariant measure of the rate of expansion. Among 
other properties, it vanishes if the horizon is non-expanding while on a dynamical horizon section it is equal to the 
rate expansion of e with respect to the unit-normalized version of the evolution vector field [39| . The examples of 
section [V] provide intuitive support for the identification of small e with a slow expansion. 

This notion of a slow area change is the key part of the entire definition. The rest of the third condition forces the 
surface to be "almost-null" with a slowly changing intrinsic metric. To simplify notation we now adapt a evolution- 
parameter-compatible scaling of the null vectors so that 



(36) 



which also ensures that any transition to equilibrium and so a fully isolated horizon will be smooth. With this scaling 
the first two conditions along with ( 23 ) imply that 



'(() 



- a ab a {l) 



T nh l a t < 



and so the rate of expansion becomes 



£ ye = —C0(n)e 

0(6») 



R\ 



(37) 



(38) 



while the "time" -derivative of the full metric is 

C vq ab = 2*$ + (C9 (n) q ab - 2Ca%> ) . (39) 



0(e 2 



The leading terms are the same as they would be for a truly null surface. Further the leading terms of the energy 
flux across the horizon 

V a T b a r b » T ab r£ b - C 2 T ab n a n b , (40) 

0(e 2 ) 0(e 4 ) 

(where r a — £ a + Cn a is the usual timelike normal) also match those of a truly null surface. 
We refer the reader to [8] for further discussion of these and other consequences of the definition. 



2. Slowly evolving horizons 

For slowly expanding horizons we imposed conditions to ensure that AH was almost null and that the two-geometry 
changed only slowly. We now move on to slowly evolving horizons where additional conditions arc imposed to also 
restrict the evolution of much of the rest of the horizon geometry. This can be motivated by the analogous development 
of the isolated horizon formalism from non-expanding horizons (restricting intrinsic geometry) to fully isolated horizons 
(also restricting extrinsic geometry). 

The two-metric q a b, extrinsic curvatures (k^j , k^t ) and normal connection uj a fully specify the geometry of the 
S v . The rest of the geometry of H is then fixed by Ky and C. From the previous section, for a slowly expanding 
horizon with an adapted evolution parameter, we have a good notion of a "time" -derivative. Thus, for these remaining 
quantities we can directly bound their rates of change with respect to V a , demanding that in each case their rate of 
change be of a lower order than the original quantity. 

Definition: Let AH be a slowly expanding section of a FOTH with a compatible scaling of the null normals. Then 
it is said to be a slowly evolving horizon (SEH) if in addition 

1. ||£yu) a || , |£yKy| and |£y#(„)| < e/R% and 

2. \C V C\ <e 3 /R H . 

Note that the inclusion the CyC is new in this article as compared to earlier definitions of slowly evolving horizons. 
It is not required for the first law however, as we shall see, it is required to prove that there is always an event horizon 
candidate in close proximity to any SEH. 



IV. KEY PROPERTIES 



In this section we review two of the key properties of slowly evolving horizons. 



A. Laws of Mechanics 



Slowly evolving horizons are intended to be the black hole analogues of quasi-static systems in thermodynamics. 
As such we would expect them to (approximately) obey the laws of black hole mechanics (a nice discussion of the 
equilibrium forms of these laws can be found in the isolated horizon chapter of this book which is also online at [2]). 

Starting with the zcroth law, wc expect that surface gravity will be approximately constant on each surface S v . 



This is the case. Applying the various defining conditions to the equation for the variation of ui a (25 1, it is immediate 
that 

K«v|| < 4- (41) 

H 

(the analogous equation for null surfaces is used to show that the surface gravity is constant for a (weakly) isolated 
horizon). Given that the surface gravity has already been required to been slowly changing up the horizon, it follows 
that over a foliation parameter range on AH that is small relative to 1/e, 

k v = k + O(e) (42) 
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for some constant n . Note however that, as for temperature in a standard quasi-static process, the surface gravity 
can accumulate larger changes over sufficiently long periods of time. 

We have already seen that FOTHs that satisfy the null energy condition arc non-decreasing in area. SEHs inherit 
this property. This is one form of the second law, however given that SEHs represent quasi-static black holes one 
might also expect them to dynamically satisfy the second law in its Clausius form. For reversible thermodynamic 
processes, entropy is defined by 

5S= 5 -^. (43) 
Switching to black hole mechanics we expect the analogous law to link area increase, energy flux and surface gravity. 



Slowly evolving horizons do obey a law of this form [8] • First linearly combining ( 23 ) and ( 24 ) one can show that 



f 



Kv 9 v = 5 v 9 {i) + B5 v e (n) + d a (Cd a C - 2Cu a ) + a^a^ + G ab V a T b + ~0 {v) B {t) , 

where the shears associated with V a and r a are defined in the obvious way. Integrating over S v eliminates the total 
derivative term and then applying slowly expanding/evolving conditions (including the zeroth law) we find that to 
order 0(e 2 ): 

e\ E inJr + T <*rn ■ ( 44 ) 



4nG k J Sv [ 8itG 

Inside the integral on the right-hand side, the terms can respectively be identified as fluxes of gravitational radiation 
and stress-energy and, as would be expected, they are equivalent to the terms that one finds for truly null surfaces 
(see for example the discussions in [9"IITB]). 

Note, though we have called this the Clausius form of the second law, it is more traditional to refer to it as the first 
law of black hole mechanics [19]. For stationary black holes and isolated horizons this first law takes the form 

SM = —— kSA + mj + $<5Q , (45) 

where the M is the black hole mass (taken as the internal energy) , J is angular momentum and Q is the charge while 
Q, and $ are their associated potentials. However, it has been pointed out [20] that while there are energy flow terms 



in (44), there is no explicit splitting of those terms into work terms versus those involving changes in the internal 



energy of the system. The notion of internal energy is key to the thermodynamic interpretations of the first law and 



as such it probably makes more sense to view the energy flow terms as more analogous to the heat flow in ( 43 ) . Hence 
we refer to it as a form of the second law. 



B. Spacetime near an SEH 



In this section we will demonstrate that (at least to leading order) there is an actual null surface close to each 
slowly evolving FOTH. This null surface is a candidate event horizon and is, in fact, the actual event horizon if the 
FOTH remains an SEH forever. Further, despite the teleological definition of event horizons the behaviour of this 
event horizon candidate is locally determined and its evolution closely mirrors that of the SEH. Full details of this 
result will appear in [21) . 



1. Event horizons 



An event horizon is the boundary of a causal black hole: a region of spacetime from which no causal signal can 
escape. Such a surface is necessarily null and, for outside observers, it is the boundary between the unobservable events 
inside the black hole and those outside that can be seen. This definition is teleological: one determines the extent (or 
existence) of a black hole by tracing all causal paths "until the end of time" and then retroactively identifying any 
black hole region. The most obvious manifestation of this type of definition is that, as shown in Figure [2] infalling 
matter curtails rather than causes the expansion of an event horizon. 

The reason for this behaviour is mathematically easy to understand. As a null surface, the event horizon is 
necessarily ruled by null geodesies. Thus, in the language of earlier sections, it is a three-surface with C — and so 
the Raychaudhuri equation holds: 



c v e w = Kt o w - (e 2 w /2 + \\a [e )\\ 2 + G ab tt b ) 



(46) 
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FIG. 2: A schematic that plots both the quasilocal and event horizons for a typical Vaidya spacetime in which a shell of dust 
(the shaded gray region) falls into a pre-existing black hole. In this figure, horizontal location measures the areal radius of 
the associated spherical shell while the direction of increasing time is roughly vertical outside the event horizon but tipping 
horizontal-and-to-the-left inside. On both sides, inward-moving null geodesies are horizontal while "outward-pointing" null 
geodesies are represented by gray dashed lines. The false event horizon is a null surface that would appear to be the event 
horizon if one is unaware that the matter shell will fall into the black hole in the future and so capture it. 



Now, given the Einstein equations and null energy condition, if we temporarily adopt an affine scaling of the null 
vectors (ki = 0), it follows that 

C e 9 w <-^e 2 w . (47) 

Thus the expansion naturally decreases with time. There is nothing profound about this statement. Consider for 
example a sphere of light expanding outwards in flat space so that R(t) = R Q + ct. Then 

"«> -*-m m 

and so the rate of expansion decreases as time (and therefore) R increases: the rate of expansion scales with the 
(inverse) radius of the horizon. Returning to general relativity a non-zero flux of gravitational energy or positive- 
energy density matter will cause the rate to decrease even more. Again this is to be expected: the gravitational 
influence of more mass inside the shell should decrease the rate of expansion. 

The apparently counterintuitive nature of event horizon evolution is then seen to be not so mysterious. Event 
horizons are families of null geodesies evolving exactly as null geodesies should. The particular set of geodesies has 
been chosen based on future boundary conditions (not escaping to infinity and not falling into the singularity) but at 
any given moment they evolve causally just like any other null congruence. 

The real mystery is then: why is it that numerical |22| and perturbative calculations |23j often show event horizons 
evolving in the apparently intuitive way, expanding in response to infalling radiation and matter? Rewriting the 
Raychaudhuri equation as 

K t e {t) -c e e (e) = \\a e f + G ab ti b + ^ef e) (49) 

it is straightforward to see that just a couple of slowly-evolving-type conditions will be sufficient to enforce a Clausius- 
type second law (a.k.a the first law of black hole mechanics). For a small 9^ one would have 0%-, <C 9(e) while one 
might also expect the rate of change of the expansion to be smaller than the expansion itself C e9^) <C (as in a 
standard derivative-expansion like those considered in the fluid-gravity correspondence |24j ) . 
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W e now consider how to express this intuition in geometrically invariant way. First note that on right-hand side of 
( 49 1 the term can be neglected if 

\olo <t(\Wtf + G ab i a e b ) , (50) 

for some e -C 1. Both sides of this inequality change in the same way with respect to scalings of the null vectors so 
this condition is geometrically meaningful and can be used to determine the scale of "slowness". Next, if we can scale 
the null vectors so that on the left-hand side 

Ct9it) ~(ib) e v>' (51) 

and K£ is of order I/Rh, then the Clausius-type law is recovered for this null surface: 

* \\a e \\ 2 +G ab e a e b . (52) 



Thus, in this limit fluxes will appear to drive the expansion of these null surfaces in the intuitive way. Further 
discussion of these slowly evolving null surfaces can be found in |25j . 



2. Event horizon candidate near an SEH 



Intuitively one might expect a truly null surface to lie close to any slowly evolving horizon. In fact, this is the 
case and that null surface is also slowly evolving in the sense discussed above. To first order this is easily seen. As 
noted earlier, the extrinsic curvature of a hypersurface is the rate of change of the induced metric under normal 
deformations. Then, the induced metric on a three-surface close to a slowly evolving horizon is 

€ T) « ?ff + P W 0) = *£> + 2(A-)*2> , (53) 
where f is the future-oriented timeline normal to AH and As measures proper time along the corresponding timelike 



geodesies. Substituting in (11) and (131 this becomes 



(Ar) 
Tab 



(2C + (A S )V2CR V ) [dv] a [dv] b + V2C(As) {Q a [dv] b + [dv] a Q b ) + (q ab + -^L + 2C*<£ ) )) , (54) 



where 

k v = k v -(C v C)/V2C^k v (55) 

(by assumption the derivative is of lower order than k\> and so may be neglected). 
Thus, the surface a proper time 

As fa ~V2C/k v (56) 

into the past of the SEH will be null to leading order. It is not hard to demonstrate that this surface is also slowly 
evolving in the sense considered in the last section and so evolves in a local way which closely mirrors that of the 
SEH 21 j. If C remains small and ultimately asymptotes to zero then this null surface also asymptotes to the SEH 
and truly is the event horizon. 

It is important to keep in mind that even in this case, the event horizon is still teleologically defined. We identify 
it based on the assumption that it will remain slowly evolving for all eternity. If, at some point in the future the 
FOTH is no longer slowly evolving, then the construction breaks down and this set of null geodesies will dive into the 
singularity as did the false horizon in Figure [2j 



V. EXAMPLES AND APPLICATIONS 



In order to build some intuition about these horizons, we now turn to a couple of concrete examples. The first 
is very simple, an SEH in a Vaidya spacetime, and it mainly serves to help to establish the range of astrophysical 
processes that we might expect to be slowly evolving. The second example is an SEH in a five-dimensonal black brane 
spacetime. This is technically much more complicated and also has the very interesting property that it is dual to a 
four-dimensional conformal fluid through the fluid-gravity duality. As such we can use it transfer lessons from slowly 
evolving horizons to fluid dynamics and vice versa. 
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A. Vaidya black holes 



We begin with a spherically symmetric example in four-dimensions: slowly evolving Vaidya black holes. This is 
taken from [26] which discusses both this example and some other, slightly more realistic, ones. As in that paper 
we are primarily interested in using the example to build physical intuition and so take take one solar mass as our 
standard unit: 

M Q = 1.9 x 10 30 kg. (57) 

Hence, with G = 6.7 x 10 -11 m 3 kg -1 s~ 2 and c = 3.0 x 10 8 ms _1 , distances and times that come out of our equations 
will be measured in units of 

R = GMq/c 2 « 1.4 x 10 3 m and (58) 
T Q = i? Q /cw 4.7 x l(T 6 s (59) 

respectively. Relative to everyday experience, our unit of mass is huge while our unit of time is very small. As such 
we might expect to find that some fairly dramatic events lie in the slowly evolving regime. We will see that this is 
the case. 



1. The spacetime 

The Vaidya solution |5] is the simplest example of a dynamical black hole spacetime. It models the collapse of null 
dust and is described by the metric 

ds 2 = - ( 1 - ^HM] dv 2 + 2d vdr + rW , (60) 



where v is an advanced time coordinate, m(v) measures the mass of the black hole on a hypersurface of constant v 
and the infalling null dust has stress-energy tensor 

T «b = (61) 

Figure [2] shows one such a spacetime in which a shell of null dust accretes onto a pre-existing black hole. 

To calculate the various quantities of interest we fix a scaling of the null vectors to surfaces of constant r and v. 



2m(v) 

1,1 U ,0,0 



and n a = [0,-1,0,0] . (62) 



Then straightforward calculations show that for these spherically symmetric two-surfaces: 



h^Ul- 2 ^ 1 ) and* (w) =-? (63) 



r v r J r 



while 

K t - -n b e a V a £ b = ^ and G ab £ a £ b = . (64) 

Overdots indicate derivatives with respect to v. Thanks to both the symmetry and the special nature of null dust, 
most other geometrical quantities (including the shear and all other components of the stress-energy tensor) vanish. 



2. The slowly evolving horizon 



Now, there is clearly a FOTH located at r = 2m(v) with characteristic length Rjj = 2m. For this surface 

C = 2m (65) 
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-CO M {R H ) 



4m 



2 ~ J (n)\ 

and it is straightforward to see that if m <C 1, then this is slowly expanding horizon. Next, 

1 



4m 



and £yC = 4m . 



(66) 



(67) 



and so, if m <C m the horizon is also slowly evolving. 

Before checking implications, we restrict to a specific example to get a feeling for the kind of situations under which 
Vaidya is slowly evolving. Consider a piecewise linear mass function : 



1 v < 

m(v) = < 1 + av < v < v a 

1 + OLV V > V 



(68) 



That is, we irradiate an initially solar mass black hole with null dust for a finite period of time. 

To interpret this physically, consider a fleet of observers who use rockets to hold themselves some distance from the 
horizon at constant areal radius r a . These observers will see a total mass of av Q fall past them during time 



T 



1 



2m(v) 



dv < av Q , 



(69) 



with the inequality coming close to saturation for very large r a . Of course, these observers can't actually see the 
horizon but since it is spacelike there aren't any suitable observers who actually live on the horizon. 

Adding in some numbers, consider a — 1/40000 and v « 2.1 x 1O 5 M (one second in standard units). Then, 
our observers would see av a w 5.3 solar masses of material sweep past them in less then a second. However, at the 
horizon a quick calculation shows that e 2 — 4a = 10~ 4 . Thus, even in this very dramatic situation the horizon would 
be slowly evolving to order e ~ 0.01. 

There are many possible objections to this example, from the fact that null dust doesn't exist to the assumption 
of spherical symmetry. However, it is in line with our expectations based on the units M Q and T Q and broadly 
comparable results can be found in [53] which includes timelike-dust examples and also Schwarzschild as perturbed by 
quite strong incoming gravitational waves. Details vary, but the general message is that black holes probably evolve 
slowly in almost all astrophysical situations, with one notable exception being black hole mergers. 



3. Mechanics 



For spherically symmetric Vaidya, the zeroth law holds trivially (and exactly) while the FOTH at r = 2m(v) 
trivially is clearly increasing in area for m > 0. Thus the only part of the mechanics that needs checking is the 
Clausius form of the second law ( 44 ) . From the preceding section 



K a dm 
8ttG ~ ~dtv 



while the right-hand side is also 



|qW|| 2 
8ttG 



dm 
dv 



(70) 



(71) 



In fact as was discussed in |26j the Vaidya spacetime is, in many ways, not the best spacetime to demonstrate that 
this law holds: thanks to the symmetries and special nature of null dust, (711 actually holds exactly whether or not 



the horizon is slowly evolving. This however is an exceptional circumstance and for all other examples (including 
Tolman-Bondi) the slowly-evolving condition is necessary. 



The event horizon candidate 



For the Vaidya spacetime one can directly search for null surfaces just outside the event horizon by trying to solve 



14 



for a solution of the form r = 2m(l + p), with p <C 1. If one assumes a hierarchy of derivatives (so that m ^> mm 3> 
m 2 m), one can iteratively solve this equation order-by-order to any desired accuracy. In particular to second order 
one finds that [24"] : 

Tnuii ~ 2m 4- 8mm + (64mm 2 + 32m 2 m) . (73) 
If this coordinate-dependent distance is converted into a proper-time measurement |25j . it matches the first order 



prediction of Eq. ( 54 1 



B. Black branes in the fluid-gravity correspondence 

A much more interesting example of near-equilibrium horizons can be found in the context of the fluid-gravity 
correspondence. This is not the place to go into a full discussion of this version of the AdS-CFT correspondence, but 
very briefly it states that a near-equilbrium asymptotically AdS black hole in A-dimensions is dual to a near-perfect 
conformal fluid in (A — l)-dimensions. From the CFT perspective, these near-equilbrium black holes correspond to 
the long-wavelength regime where the field theory reduces to fluid mechanics. That said, from a purely gravitational 
perspective, there is no need to reference the full correspondence. A direct examination of the long-wavelength pertur- 
bations of AdS black branes demonstrates that the gravitational equations of motion governing those perturbations 
arc equivalent to the Navier-Stokes equations governing a near-perfect conformal fluid |24j . 

Thermodynamically the temperature of this fluid should correspond to that of the black brane while the entropy 
flow should map onto the expansion of its horizon, and indeed this is the case to leading order. However, at higher 
order things become more interesting and several ambiguities arise in the exact meaning of the phrase "the entropy 
flow should map onto the expansion of its horizon" . 

In this section we briefly review the role of slowly evolving horizons in the fluid-gravity correspondence. Further 
details can be found in |27H29) . Up until now we have restricted our attention to four-dimensional spacetimes with 
spatially compact black holes. Switching to five-dimensions for this example, the formalism that we have developed 
remains essentially unchanged. The only changes are that we relax the compactness condition (with the scale now 
being set by the surface gravity rather than the areal radius) and have to replace some factors of 1/2 with 1/3 in 
equations such as (j6j, ([23]) -(26} and ((27 1. 



1. The spacetime 

Starting with a five-dimensional AdS black-brane metric, the fluid gravity correspondence may be derived in the 
following way |24| . This solution is static and on inducing a compatible time- foliation, each instant may be thought 
of as a stack of three-dimensional intrinsically flat planes that stretch from the spacetime singularity to infinity with 
a similarly planar event horizon in between. In Eddington-Finkelstein-like form the metric is 

ds 2 = -r 2 ( 1 - —^—j) dv 2 + 2dvdr + r 2 (dx 2 + dy 2 + dz 2 ) , (74) 



(rb o y 

where for computational convenience the cosmological constant has been scaled to take the value —6 and the event 
horizon at r = l/b a has temperature 

T= ^h- (75) 

An alternate form of this metric arises if one performs a Lorentz boost within the planes so that we switch to a 
coordinate frame which moves with four-velocity u" relative to the static system: 

ds 2 = -r 2 ( 1 l — ) u u v dxi*dx v + 2u°dx f *dr + r 2 {n„ v + u°u° v ) dx^dx" . (76) 

V ( rb °) J 

In this case the greek indices run over {v,x,y,z} and rj^ is the standard four-dimensional Minkowski metric. 

The fluid-gravity correspondence comes from considering long-wavelength perturbations of this boosted black brane 
metric: one replaces the constant b oy u" and rjau with b, u a and g^ v which agree with their unperturbed counterparts 
at leading order but beyond that have a coordinate dependence. In the planar directions the scale of variation of these 
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perturbations must be much longer than the length scale set by b (they are long- wavelength perturbations). Thus, 
for example, 



WdaU^W < 



(77) 



Further, one assumes that this pattern continues for the entire hierarchy of derivatives in the planes. Continuing with 
the u 13 example: 

\\u s \\ > b\\d a u s \\ > b 2 \\d a d p u 5 \\ > b 3 \\d a d d jU s \\ . . . (78) 

It is then possible to separate and solve the Einstein equations for these perturbations order-by-order. As might be 
expected, though straightforward in principle, these computations rapidly become very intricate [2~31 I33j . 

For those interested in horizons these approximate solutions are of interest in their own right. At second order 
in this perturbation theory we can track the evolution of both event and apparent horizons and they are distinct 
from the unperturbed results. By contrast for Schwarzschild spacetimes current state-of-the-art perturbation theory 
cannot distinguish these horizons and places them both at r rs 2m; corrections to the position are at a higher order 
than can be tracked [30 32 . However, before considering those horizons in more detail we complete our discussion of 
the fluid-gravity correspondence as it turns out that the conformal fluid interpretation eases the calculations on the 
gravity side. 

The perturbed solution is entirely described by the quantities {b, u^, g^} with all tensor-indices being in the 
cotangent space to the surfaces of constant r. Each of these surfaces has an intrinsic geometry defined by g^ v (it 
is weakly perturbed from Minkowski) and on each such surface one can construct a four-dimensional stress-energy 
tensor: 

= eu^u v +ph^ , (79) 

where 

h^u = g^u + u^Uy (80) 

is the three-metric transverse to the and the pressure and energy density are defined by b and obey the equation 
of state: 

£ = 3 ^16^- < 81 > 

This may be regarded as a formal definition (which will be justified by later developments), however it may also be 
derived systematically as a Brown- York boundary stress-energy tensor [35]. In any case, for the unperturbed metric 
this is the stress-energy of a conformally invariant perfect fluid. However on implementing the perturbations it is 
only near-perfect with the corrections depending on gradients of {&, u 11 , g^ v } (and corresponding to non-perfect-fluid 
corrections such as viscosity). 

Implementing the Einstein equations turns out to be equivalent to requiring conservation of the stress-energy tensor 
on each of these surfaces: 

Vf/T^" = , (82) 

where T>^ is the conformally invariant covariant derivative compatible with g^ u (see the original sources or |29) for 
more details). This can be observed directly from computation with the components, but it may also be seen as a 
consequence of the fact that the momentum constraint must hold on each surface (again see [35 ) . 

The dual fluid is generally taken to live on the surface at r — > oo. The conformal invariance can be used to rescale 
away diverging quantities in this limit and so guarantee that the theory remains well-defined. 



2. The slowly evolving horizon 



For the Vaidya example, finding a FOTH was straightforward; guided by the spherical symmetry, we simply 
examined the surfaces of constant r (and ignored all other candidates). Here things are a bit more complicated. The 
geometric symmetry of the problem is broken and so one cannot use it to select and search for a preferred FOTH. 
However, there is still the conformal symmetry and this turns out to be sufficient to select a preferred FOTH (to 
second order). Again we will be brief, leaving the reader to consult [2H UHl 133 GH] for more details. 
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The conformal symmetry of the fluid means that solutions are invariant under rescalings of the form 

g^v^e-^g^, u^^e^u", b^e'H and r -> e*r . (83) 

Physically meaningful quantities should also respect this symmetry and we take this to include the location and 
evolution vector field of the FOTH. This is a strong restriction as there are only a limited number of conformally 
invariant quantities available at each order of the gradient expansion. For example, at first order, there is only the 
Weyl-invariant shear of u^: 

a^ d = V ( ,u v) , (84) 

where the parentheses indicate the usual symmetrization. At second order in gradients there are ten quantities (three 
scalars, two vectors, and five transverse traceless tensors). We require that corrections to both the location of the 
horizon and its evolution vector field V M be constructive out of these quantities. 

After some fairly extensive algebra one can show that to second order there is a unique conformally invariant FOTH 

at 

™- i + i (*."*- (85) 

where h^ H is a non-trivial numerical factor, while the conformally invariant factors are built out of the conformally 
invariant shear and vorticity of as well as the similarly conformally invariant curvature scalars: S\ = b 2 a uu <j^ , 
S% = b 2 uj llv Lo^ 1 ' and S3 — b 2 lZ. Note that the corrections are all at second order in gradients. At first order the 
horizon is still at r = 1/6. This FOTH is easily seen to be slowly evolving with characteristic length scale Rr = l/b. 

The exact implications of this uniqueness are unclear. It has been conjectured |29j that the uniqueness might extend 
to higher order corrections and possibly even imply that in this case the conformal invariance selects a unique FOTH 
just as in the Vaidya case the spherical symmetry selected a preferred horizon. If true this would be a fascinating 
result, however it is also possible that the uniqueness to second order is simply a consequence of the fact that the 
horizon is slowly evolving. The resolution of this question will be left to future investigations. 



3. Mechanics: horizon and fluid 



On the gravity side it is easy to check that both the (approximate) zeroth law holds («y « l/b) as does the Clausius 
form of the second law. However for these solutions there is an additional interest as we can compare the mechanics 
of the near-equilibrium horizon to the thermodynamics of the similarly near-equilbrium fluid. By the fluid-gravity 
correspondence we would expect them to match. 

Before we can do this, we first need to understand the thermodynamics of the fluid. It is standard to identify the 
fluid temperature with that of the horizon so that to first order 

1 "fluid ~ t • (86) 



Then the equation of state (81) defines energy density and pressure as functions of temperature. 

The next step is to identify an entropy current J M : the divergence of such a current V^J^ determines the rate of 
entropy production in the fluid. For equilibrium flows this takes the form 

J& = su» (87) 

where s is the (equilibrium) entropy density. However, away from equilibrium the correct form of the entropy current 
is not so clear and it is usually argued that one should allow for correction terms [MJ [STJ I2H ISS] ■ In our situation 
we must consider all conformally invariant causal flows with non-negative divergence. As for the horizons, requiring 
that the flow be conformally invariant significantly restricts the freedom of construction, but in this case is does not 
completely remove it: there are allowed second order corrections to J£ q and these correspond to different values of 
entropy production. 



The rate of expansion of the slowly evolving FOTH defined by ( 85 ) corresponds to one of these values of entropy 
production. Given that FOTHs are non-unique (recall that they are not rigid and may be deformed), it is very 
tempting to then conjecture that the other J M correspond to other (nearby) FOTHs. However, this cannot be the 
case: even though the conformal symmetry does not select a unique entropy current, it did select a preferred slowly 
evolving horizon. Thus any gravitational interpretation of this uncertainty must be sought elsewhere. 
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4- The event horizon candidate 

Our FOTH is slowly evolving, and so we would expect to find a (similarly slowly evolving) event horizon candidate 
in close proximity. This is the case [213 anc ^ so this candidate provides another potential gravitational dual to the 
entropy current. Indeed it is straightforward to check that the rate of expansion of the event horizon candidate is 
matched by the divergence of another member of the family of allowed entropy currents. 

5. Implications and complications 

From the horizon perspective, these black brane spacetimes with their (measurably) dynamical horizons are nice ex- 
amples in their own right providing testing grounds for ideas about slowly evolving horizons. The fluid gravity duality 
adds another level of interest to these examples. In the first place we can see that dual to the near-equilibrium horizons 
are the similarly near-equilibrium fluids and that the mechanics of the horizon is similarly dual to a thermodynamics 
of the fluid. 

Note the phrase "a thermodynamics" rather than "the thermodynamics" . The entropy flow in the fluid is not 
uniquely defined and so therefore neither is its thermodynamics. Both the FOTH and event-horizon candidates 
have their own (closely-related but distinct) thermodynamic duals. However, these are just two of the potential 
thermodynamics and so it is natural to consider duals to the other entropy flows. There are at least two possible ways 
to account for this freedom in a geometric way: 

1. The conformal fluid lives on the AdS boundary of the spacetime and should be dual to the entire spacetime, not 
just the horizon. It can be argued that the brane dynamics dominate the spacetime, however even in that case one 
must still define a map between the horizon and the boundary. Given the coordinate system, it is computationally 
easiest to just match up points with matching coordinates. However this mapping is not geometrically 
preferred and there are other acceptable ways to set it up. It turns out that the freedom in constructing the 
mapping closely matches that in defining the entropy flow[24 . So, the first possible geometric explanation for 
the entropy flow freedom is that it is a manifestation of the horizon-boundary mapping uncertainty. Either the 
FOTH or the event horizon may be viewed as dual to each of the entropy flows. 

2. The second possibility is to take the freedom on the fluid side as signalling that a similar freedom should exist 
in defining entropy and entropy flows in the bulk. References [271 128] suggest one way of doing this in which the 
SEH and event horizon candidate are taken to be just two examples of a family of horizon-like structures any 
one of which can be associated with a mechanics (and each of which is dual to a fluid flow). These surfaces can 
have any signature but each signals the existence of a nearby casual boundary and is (at least approximately) 
null and slowly expanding with Oig\ ps 0. In any approach to equilibrium they all are required to have the same 
isolated horizon limit. In particular the timelike ones are similar to those surfaces considered in the membrane 
paradigm [37 . 

Of course these two possibilities are not mutually exclusive. If the conformal fluid reflects the true physics of the system, 
then the first resolution means that the horizon-candidates proposed by the second resolution are indistinguishable 
since both have the same range of freedom in their definition. Thus each candidate would furnish an equally acceptable 
dynamical black hole mechanics. These would differ in detail but remain closely related and in particular all reduce 
to the same equilibrium limit. 

A definitive resolution of these issues will have to wait for future work. However it may be worth pointing out that 
uncertainties of this type would not be too surprising. While there are several well-defined measure of total energy in 
general relativity pQ , its localization is famously ill-defined and the study of potential quasi- local definitions of energy 
remains an area of active research [38j . Thermodynamics and black hole mechanics are intimately tied to energy 
flows and so if there are uncertainties in amounts of energy transferred it would be consistent if there were similar 
uncertainties in entropy. 

VI. SUMMARY 

The first law of black hole mechanics as formulated for exact solutions or isolated horizons is a phase space relation: 
it considers variations through the phase space of equilibrium states rather than a dynamical evolution between those 
states. The original motivation for the investigation of slowly evolving horizons was to better understand the near- 
equilibrium regime for dynamical trapping horizons and so reinterpret this first law in a dynamical way [7] . We have 
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seen that this reformulation is possible and the resulting dynamical version arises as a particular limit of the equations 
governing the geometric deformation of the surfaces. 

For null surfaces the deformation equations reduce to the Raychaudhuri equation and so also govern the evolution 
of event horizons, though the identification of any particular null surface as the event horizon remains a teleological 
procedure. There is also a slowly evolving limit for null surfaces and we saw that there is a surface of this type in 
close proximity to each slowly evolving horizon. If the FOTH remains near-equilibrium at all times in the future, that 
null surface will be the event horizon. Its physical properties and geometry are very similar, though not identical, to 
those of the slowly evolving horizon. 

A wide variety of astrophysical processes are probably in the near-equilibrium regime. One can quickly see this by 
noting that the characteristic time scale for black hole would be the horizon-crossing time. For a solar mass black 
hole this is on the order of 10~ 6 seconds. Thus we would expect anything that happens on a time scale much greater 
than this to be slowly evolving. Particularly dramatic examples were demonstrated in the Vaidya spacetime. 

Finally we considered the interesting case of five-dimensional black branes and the fluid-gravity duality. In that 
case the perturbatively constructed spacetime gives rise to a slowly evolving horizon in the bulk and a corresponding 
near-equilibrium fluid on the boundary, with the mechanics of the horizon closely matching the thermodynamics of 
the fluid. However, the thermodynamics is not uniquely defined and can be adjusted to match both the SEH and 
its attendant event horizon candidate. It can be matched with the evolution of other nearly null, nearly isolated 
"horizons". This behaviour is reminiscent of the classical membrane paradigm which noted that much black hole 
physics can be modelled by considering the interaction of near-horizon timelike surfaces with their environment. 

Acknowledgements 

This work was supported by the Natural Sciences and Engineering Research Council of Canada. 



[1] R.M. Wald, General relativity (University of Chicago Press, 1984). 

[2] J. Engle and T. Liko, "Isolated horizons in classical and quantum gravity", arXiv:1112.4412. 

[3] A. Ashtekar and B. Krishnan, "Isolated and dynamical horizons and their applications, Living Rev. Rel. 7, 10 (2004). 
[4] A. Ashtekar, C. Beetle, O. Dreyer, S. Fairhurst, B. Krishnan, J. Lewandowski and J. Wisniewski, "Generic isolated horizons 

and their applications", Phys. Rev. Lett. 85 3564 (2000). 
[5] E. Gourgoulhon and J. L. Jaramillo, "A 3+1 perspective on null hypersurfaces and isolated horizons", Physics Reports, 

423 159 (2006). 

[6] I. Booth, "Black Hole Boundaries" Can. J. Phys. 83 1073 (2005). 

[7] I. Booth and S. Fairhurst, "The first law for slowly evolving horizons", Phys. Rev. Lett. 92, 011102 (2004). 

[8] I. Booth and S. Fairhurst, "Isolated, slowly evolving, and dynamical trapping horizons: Geometry and mechanics from 

surface deformations", Phys.Rev, D75 084019 (2007). 
[9] E. Poisson, A Relativist's Toolkit: The Mathematics of Black- Hole Mechanics (Cambridge University Press, 2004). 
[10] S.A. Hayward, "General laws of black hole dynamics", Phys.Rev. D49 6467 (1994). 

[11] S.W. Hawking and G.F.R. Ellis, The large scale structure of space-time, (Cambridge University Press, 1973). 
[12] A. Ashtekar, S. Fairhurst and B. Krishnan, Phys.Rev. D62 104025 (2000); A. Ashtekar, C. Beetle and J. Lewandowski, 
Phys.Rev. D64 044016 (2001). 

[13] A. Ashtekar and B. Krishnan, "Dynamical horizons: energy, angular momentum, fluxes and balance laws", Phys. Rev. Lett. 

89 261101 (2002); "Dynamical horizons and their properties", Phys. Rev. D68 104030 (2003). 
[14] A. Ashtekar and G. J. Galloway, Adv. Theor. Math. Phys. 9, 1 (2005) [gr-qc/0503109]. 
[15] I. Bendov, "The Penrose inequality and apparent horizons", Phys.Rev. D70 (2004) 124031. 

[16] I. Booth, L. Brits, J. Gonzalez, and C. Van Den Broeck, "Marginally trapped tubes and dynamical horizons" 

Class. Quant. Grav. 22, 4515 (2005). 
[17] I. Booth, "Two physical characteristics of numerical apparent horizons", arXiv:0709.0934 (2007). 
[18] S. W. Hawking and J. B. Hartle, Commun. Math. Phys. 27 283 (1972). 

[19] J. M. Bardeen, B. Carter, S. W. Hawking, "The Four laws of black hole mechanics", Commun. Math. Phys. 31 161 (1973). 
[20] E. Gourgoulhon and J. L. Jaramillo, "Area evolution, bulk viscosity and entropy principles for dynamical horizons," Phys. 

Rev. D 74, 087502 (2006). 
[21] I. Booth, "Spacetime in the vicinity of a slowly evolving horizon" (in preparation). 

[22] P. Anninos et. al, Phys. Rev. Lett. 71 2851 (1993); P. Anninos et. al., Phys. Rev. Lett. 74, 630 (1995); Baiotti et.al, 

Phys. Rev. D 71 024035 (2005); Schnetter et. al. Phys. Rev. D 74 024028 (2006) 
[23] E. Poisson, Phys. Rev. D 70 084044 (2004); E. Poisson, Phys. Rev. Lett. 94 161103 (2005); E. Poisson and I. Vlasov, 

Phys. Rev. D 81 024029 (2010). 



19 



[24] S. Bhattacharyya, V. E. Hubeny, S. Minwalla and M. Rangamani, " Nonlinear Fluid Dynamics from Gravity" JHEP 0802, 
045 (2008); S. Bhattacharyya et al, "Local Fluid Dynamical Entropy from Gravity" JHEP 0806, 055 (2008). 

[25] I. Booth and J. Martin, "Proximity of black hole horizons: Lessons from Vaidya spacetime", Phys.Rev. D82 124046 (2010). 

[26] W. Kavanagh and I. Booth, "Spacetimes containing slowly evolving horizons", Phys.Rev. D74 044027 (2006). 

[27] I. Booth, M. Heller and M. Spalinski, "Black brane entropy and hydrodynamics: The boost-invariant case", Phys. Rev. 
D80 126013 (2009). 

[28] I. Booth, M. Heller and M. Spalinski, "Black brane entropy and hydrodynamics", Phys. Rev. D83 061901 (2011). 
[29] I. Booth, M. Heller, G. Plewa and M. Spalinski, "On the apparent horizon in fluid-gravity duality", Phys. Rev. D83 106005 
(2011). 

[30] E. Poisson, "Retarded coordinates based at a world line, and the motion of a small black hole in an external universe", 
Phys. Rev. D69, 084007 (2004); "Absorption of mass and angular momentum by a black hole: Time-domain formalisms 
for gravitational perturbations, and the small-hole/slow-motion approximation", Phys. Rev. D70, 084044 (2004); "Metric 
of a tidally distorted, nonrotating hole", Phys. Rev. Lett. 94, 161103 (2005). 

[31] E. Poisson and I. Vlasov, Phys. Rev. D 81, 024029 (2010) 

[32] I. Vega, E. Poisson and R. Massey, "Intrinsic and extrinsic geometries of a tidally deformed black hole," Class. Quant. 
Grav. 28, 175006 (2011) 

[33] S. Bhattacharyya et. al, "Conformal nonlinear fluid dynamics from gravity in arbitrary dimensions", J. High Energy Phys. 
12 116 (2008). 

[34] R. Loganayagam, "Entropy Current in Conformal Hydrodynamics", JHEP 0805, 087 (2008). 

[35] J. D. Brown and J. W. York, Jr., "Quasilocal energy and conserved charges derived from the gravitational action," Phys. 
Rev. D 47, 1407 (1993). 

[36] P. Romatschke, "Relativistic Viscous Fluid Dynamics and Non-Equilibrium Entropy," Class. Quant. Grav. 27, 025006 
(2010) 

[37] K. S. Thorne, R. H. Price and D. A. Macdonald, Black holes: the membrane paradigm, (Yale University Press, 1986). 
[38] L. B. Szabados, "Quasi-Local Energy-Momentum and Angular Momentum in GR: A Review Article," Living Rev. Rel. 7, 
4 (2004). 

[39] One might be concerned as to whether this quantity vanishes in the limit C — > (in which case V a is no longer defined). 
It does. For an explicit demonstration see the discussion in [17] . 



